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Abstract: We study spin 1/2 isoscalar and isovector candidates in both parity 
channels for the recently discovered B + (1540) pentaquark particle in quenched lattice 
QCD. Our analysis takes into account all possible uncertainties, such as statistical, 
finite size and quenching errors when performing the chiral and continuum extrap- 
olations and we have indications that our signal is separated from scattering states. 
The lowest mass that we find in the I p = 0~ channel is in complete agreement with 
the experimental value of the G + mass. On the other hand, the lowest mass state in 
the opposite parity I p = + channel is much higher. Our findings suggests that the 
parity of the + is negative. 
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1. Introduction 

The possible existence of exotic hadrons has long been put forward in several con- 
texts, but the experimental discovery of the first such particle had to wait until early 
this year 0-0. For a long time exotic hadron states containing more than three 
quarks were not considered to be of great importance on account of their presumed 
large decay widths. (For a recent discussion of this issue see e.g. f5[.) The experimen- 
tal discovery of the 9 + (1540) particle changed this situation dramatically. Indeed, 
the tightest experimental upper bound so far on the width of the 9 + is around 
10 MeV 0], and based on the reanalysis of older experimental data, more stringent 
constraints on the width are also suggested || . This remarkably narrow width would 
explain why the G + has not been seen before. 

Since the G + decays into a neutron and a K + , its strangeness has to be +1, 
the third component of its isospin is 0, and its minimal quark content is dduus. 
From the lack of a signal in the 1% = 1 channel the SAPHIR collaboration concluded 
that the 9 + should be an isospin singlet state ||. Its spin and parity cannot be 
unambiguously pinned down based on currently available experimental data. 

The experimental discovery of the + pentaquark triggered a flurry of theoretical 
speculations J7| about its possible structure, yet unmeasured quantum numbers and 
on the possibility of the existence of other exotic hadrons. Chiral soliton models M, 
as well as correlated quark models ||, [10], |Tl[ suggest that the 9 + should be a spin 



1/2 positive parity isosinglet state. In contrast, uncorrelated quark models predict 
negative parity. QCD sum rules have also been used to study the © + [|12]] . 

In the present work we use lattice QCD to study the properties of some pen- 
taquark states that are likely candidates to be identified with the + . This is the only 
available tool to extract low energy hadronic properties starting from first principles, 
i.e. QCD. In our exploratory study we choose to work in the quenched approxima- 
tion, which is known to be successful in reproducing mass ratios of stable hadrons 
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Working at three different values of the lattice spacing, a = 0.17, 0.12 and 0.09 fm 
enables us to do a continuum extrapolation. Since using light enough quarks that 
reproduce the physical pion mass would be prohibitively costly, for each lattice spac- 
ing we simulate at four different, somewhat larger pion masses and extrapolate to 
the physical limit. This is done in two different ways, by keeping the s-quark mass 
fixed to reproduce the physical kaon mass in the chiral limit and also by changing 
the s-quark mass in a way to keep the nucleon to kaon mass ratio fixed as the light 
quark mass is varied. In order to have finite size effects under control we always keep 
the linear physical size of the spatial box fixed at L « 2 fm. In addition, on the two 
coarsest lattices, we also perform a finite volume analysis by considering different 
volumes. 

Since the exact structure of the G + is yet unknown, there are quite a few differ- 
ent possibilities to construct a baryon with the given quark content dduus. In our 
computation we use two particular spin 1/2 operators, one with isospin 1 = and 
and another with 1=1. The analysis is carried out in both parity channels. 

The main result of the present work is that the lowest mass that we find in the 
I p = - channel is in complete agreement with the experimental value of the G + 
mass. On the other hand, the lowest mass in the opposite parity I p = + channel 
turns out to be much larger, thus we suggest negative parity for the + particle. 

The lightest state in the 1 = 1 channel also has negative parity and comes out 
to be about 15% heavier, but still within one a of the G + mass. The SAPHIR 
Collaboration, however, have already ruled out the + to be a member of an / > 
multiplet ||, therefore the triplet state that we found might be a genuinely new state 
not seen before, albeit with a potentially much larger width than that of the + . 

The rest of the paper is organized as follows. In Section ^ we give the details 
of our simulations and our data analysis. In Section [| we present our results and 
conclusions. Those who are not interested in lattice technicalities can directly skip 
to Section 0. 



2. Simulation details and data analysis 

In this section we discuss some technical details of our lattice simulations. We use the 
Wilson gauge action and the Wilson fermion action. The only "non-standard" piece 
of the calculation is the choice of the pentaquark operator. Since there is general 
agreement in the literature that the G + is very likely to have spin 1/2, we restricted 
ourselves to this choice. Nevertheless, it will be interesting to look at higher spin 
states too. After fixing the quark content of the particle to be dduus and its spin 
to 1/2, there is still considerable freedom in choosing the interpolating operator. 
Even though the + is experimentally suggested to be an isosinglet, we computed 
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correlators both in the I 



Vo/i 



and 1 = 1 channels using the interpolating fields 
e ahc [u T a C lb d b ]{u e s e i lb d c T ( u ^d)}. (2.1) 



Here C is the charge conjugation operator, only colour indices are explicitly written 
out and the 1 = case corresponds to the minus sign. There are many more possi- 
bilities even in the 1=0,1 channels. In principle testing other interpolating operators 
and checking which one has the best overlap with the G + state should provide infor- 
mation on the wave function of the particle. In this exploratory study, however, we 
do not pursue this direction any further. 

We did not study the 1 = 2 channel either, mainly because it seems to be com- 
putationally more expensive. At first sight this is hard to understand since normally 
in lattice simulations one is used to quark propagator computations being the most 
expensive part. Given the propagators, calculating the different contractions en- 
tering hadron correlators takes only a tiny fraction of the time. With these more 
complicated pentaquark hadrons, however, that is no longer the case. Already in the 
I = 0, 1 channels the time for computing the contractions is of the same order of 
magnitude as the fermion matrix inversion. The reason for that is mainly that the 
number of terms occurring in the colour and Dirac index sums grows exponentially 
with the number of constituent particles. 

The + can decay into a kaon 
and a nucleon even in the quenched 
approximation. The experimentally 
very narrow pentaquark state is just 
slightly above the N + K threshold, 
thus one expects that the + is vis- 
ible in correlators of our well localised 
operators. Note that in principle scat- 
tering states can mix in the spec- 
trum, though they have a charac- 
teristic volume dependence. (This 
volume dependence can be used to 
calculate decay rates c.f. fT?!). On 
the two coarsest lattices we searched 

for any possible volume dependence of the G + by considering different spatial lat- 
tice sizes ranging from L = 1.3 to 2.7 fm. We did not see the expected volume 
dependence characteristic to the S-wave scattering state. (Note, however that this 
signature would appear only in the 1 = 1 channel, since in the 1 = channel the 
scattering length is too small.) In addition, on our largest lattice at the largest quark 
mass, based on 300 configurations, we computed the 2x2 correlation matrix of the 
rj operator of Eq. ( |2.1|) and a similar operator with its colour indices contracted as 
in the nucleon and the kaon. The latter is expected to couple more strongly to the 
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0.093 


1.86 


94 



Table 1: The seven gauge field ensembles used 
for the study of the pentaquark. The lattice 
spacing has been set with the Sommer parame- 
ter ro = 0.50 fm. 
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scattering state. After analysing the cross correlator matrix we obtained for the 
smallest energy state E /(mN + rrix) = 0.994(18) (which we interpret as the S-wave 
scattering state) and for the state above Ei/(rriN + txik) = 1.074(20) (which is the 
@ + candidate). Note that the volume of our spatial box implies that the next scat- 
tering state with zero total momentum should not occur below E /(rriN + mK) = 1.21. 
Thus, in particular, we do not expect to be faced with the unlucky scenario that in 
our analysis one of the scattering states is accidentally lifted to the + mass. 

Given our choice of interpolating operators for the pentaquark, we now discuss 
other details of the simulations. We generated seven different ensembles of quenched 
gauge field configurations with the Wilson plaquette action. Quark propagators 
were calculated using Wilson fermions. The parameters of our gauge ensembles are 
summarized in Table [1|. Ensembles with three different values of the coupling were 
chosen to have roughly the same spatial volume in order to avoid any finite size 
contamination. 

The last ingredients of our computation are the choice of quark masses and the 
way we fix the two independent parameters, the lattice spacing and the s-quark 
mass. Due to the extensive literature on quenched spectroscopy one can safely plan 
the necessary parameter sets in order to carry out acceptable chiral and continuum 
extrapolations. We used light quark masses in the range m n /m p = 0.5 — 0.7. In 
quenched QCD there is a small inherent ambiguity in how those parameters can 
be fixed. Given the fact that the + is intimately connected to the kaon, it seems 
reasonable to use the physical mass of the kaon to set the s-quark mass. We adopted 
two different schemes in approaching the chiral limit. In both cases the light quarks 
were taken degenerate. In scheme (A) we kept the s-quark mass fixed to reproduce 
the physical kaon mass in the chiral limit. In scheme (B) we changed the mass of the 
s-quark along with that of the light quarks so as to keep the nucleon to kaon mass 
ratio constant. 

One of the main questions we set out to address in the present study is to 
determine the parity of the + particle. Simple baryon interpolating operators in 
general are expected to couple to states of both parities |I8[ fjl|. With the parity 
assignment of our operators the correlators are of the form 1 

(B B t ) = (1 - 7o ) [C+e-' m + + C.be-^ m -] + (1 + 7o ) [C + be~^ m + + C_e" im -] , 

(2.2) 

where T is the size of the lattice in the time direction, b = ±1 corresponds to 
periodic/antiperiodic boundary conditions in the time direction, and m± are the 
masses of the ± parity states. To be able to project out the two parities separately, 

1 In the first version of this paper the signs of the terms proportional to 70 were opposite. We 
thank S. Sasaki (2(J and K.F. Liu (2l] for pointing out this apparent error. The only consequence of 
correcting this error is the flip of the parity assignment. All the other results, including correlators 
and masses remained unchanged. 
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we always computed both components of the correlators, i.e. the one proportional to 
1 and 7o- In addition, at one quark mass on each ensemble we also computed the 
correlators with antiperiodic boundary condition, which enabled us to determine the 
two parity partners in another way. 

We performed two-parameter correlated fits to the 1 and 70 components of the 
(BoBt) pentaquark correlators with the standard cosh and sinh form respectively, 
assuming that at large enough separation the lightest state dominates the correlators. 
Since our data sample was rather limited we followed Ref. [2^j and used correlated 
fitting with smeared smallest eigenvalues of the correlation matrix. We could always 
find fit ranges with acceptable x 2 ~ 1 accompanied by a plateau in the effective and 
in the fitted masses. The fits to the 1 and 70 components always gave compatible 
results. The interpolating operators, however, could be better chosen, since the 
mass plateaus were reached at relatively large time separations compared to e.g. the 
nucleon correlator. In Fig. [I] we show a typical set of fits in the 1 = channel. 
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Figure 1: A typical fit for a correlator in the 1 = channel at j3 = 5.85 at our lightest u 
quark. The mass given by the fit along with the corresponding \ 2 P er degree of freedom is 
shown versus the starting point of the fit range. 

Since one of our most important goals is to determine the parity, let us discuss 
this issue in some detail. The interpolating operators we used are expected to couple 
to states of both parities in the given channel. In order to determine the parity of 
the lowest state as well as the mass of its heavier parity partner we performed single 
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exponential fits to both the 1 + 70 and the 1 — 70 components of the correlators. 
As can be seen from Eq. (|2.2|) these projections should fall off with two different 
exponents for < t <C T/2 and for T/2 <C t < T, corresponding to the two different 
parity states. As expected, we could always confirm that one of the exponents was 
compatible with our previously found lowest mass obtained by the cosh and sinh fits 
to the 1 and 70 components, respectively. Both for the 1 = and the 1=1 channels, 
the negative parity state turned out to be the lighter one of the two parities. 

The above correlators were calculated with periodic quark boundary conditions. 
An independent method to determine the parity is based on the combination of the 
correlation functions obtained with periodic and antiperiodic boundary conditions. 
We carried out this analysis at each lattice spacing for the largest quark masses. As 



can be seen from Eq. (|2.2|) , the sum (difference) of the 1 component of the peri- 
odic and the 70 component of the antiperiodic correlator projects onto the negative 
(positive) parity state, respectively. Single mass fits to these combinations yielded 
masses in complete agreement with our mass and parity determinations using peri- 
odic boundary conditions. 

A consistent picture emerged for the parity of the pentaquark states at all lat- 
tice spacings and quark masses. The negative parity mass was always by many 
standard deviations below the positive parity mass. This result was obtained from 
both boundary conditions combining the 1 and 70 channels and also from combining 
different boundary conditions. Thus, we conclude that both for the 1 = and the 
1 = 1 channels, the negative parity state was determined to be the lighter one of the 
two parities. 

The masses in the I p = O 1 * 1 , l 1 * 1 channels were then extrapolated to the chiral 
limit. As we already mentioned, for the chiral extrapolations we used two different 
sets of quark masses and correspondingly two different methods. 

(A) We set the scale by the Sommer parameter ro = 0.5 fm and we held the s 
quark mass fixed to reproduce the physical kaon mass mx = 494 MeV. This 
approach resulted in strange quark hopping parameters of k s =0. 1621,0. 1574 
and 0.1544 for /3=5. 70,5.85 and 6.00, respectively. At all three lattice spacings 
four different light quark masses were used in order to extrapolate into the 
chiral limit. 

(B) A somewhat different method to approach the chiral and continuum limit can 
be imagined as follows. For all the three lattice spacings we used four sets of 
strange and light hopping parameters. The light quark masses were used to 
extrapolate into the chiral limit, whereas we chose the strange quark mass to 
reproduce the experimental nucleon to kaon mass ratio. 

We followed both methods all the way through, first performing the chiral extrapo- 
lation at fixed j3 then carrying out the continuum extrapolation. The two methods 
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(A and B) should give identical final results, though the difficulties connected with 
the chiral and continuum extrapolations might be quite different. Luckily enough, 
the extrapolations for the relevant mass ratios were equally easy for both techniques. 
The obtained results are in complete agreement. Since we used the same gauge con- 
figuration ensembles for the two methods, the results are strongly correlated and 
cannot be combined. In what follows, we focus on the somewhat more conventional 
method (A). 
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Figure 2: A typical chiral extrapolation of hadron masses on our finest lattice at (3 = 6.0. 
The particles are (from top to bottom) the 1=1,0 negative parity pentaquark states, the 
nucleon and the kaon. 

The chiral extrapolations for hadron masses can be performed with the fit func- 
tions 

m 2 H = a + bm 2 or m 2 H = a + bra 2 + cm\. (2-3) 

The difference between the two extrapolations gives some information about sys- 
tematic uncertainties in the extrapolated quantities. Note, however, that our quark 
masses are quite small and we have only four different quark masses at each (3. Thus 
with our statistics the quadratic fit turned out to be appropriate. Comparing our 



spectrum with similar results in the literature (c.f. [14|) we concluded that the un- 
certainties due to the chiral logarithms in the physical limit were subdominant at 
our present statistics. A typical chiral extrapolation for the kaon, the nucleon and 
the negative parity I = 0, 1 pentaquark states is shown in Fig. 0. Note, that our 
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chirally extrapolated masses are merely illustrations. Since quenched spectroscopy 
is quite reliable for mass ratios of stable particles, it is physically even more moti- 
vated to extrapolate mass ratios instead of masses (or equivalently one might use the 
squared ratios suggested by the above mass fit formula). As we will see the quark 
mass dependence of these ratios can be weaker than that of the individual hadron 
masses. 

The pentaquark states have the same quark content as the nucleon and the 
kaon. Therefore, a particularly attractive dimensionless quantity to look at is the 
mass ratio of a pentaquark state to the nucleon plus kaon mass. We also did chiral 
extrapolations of this mass ratio squared, and an example of that for both isospins 
and parities is shown in Fig. H. 
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Figure 3: Chiral extrapolation of the squared mass ratios m| q / (m^r + m^) 2 on our finest 
lattice at f3 = 6.00. 



Finally, we performed a continuum extrapolation for the chirally extrapolated 
quantities. In Fig. [| we demonstrate the continuum extrapolation of the pentaquark 
to nucleon plus kaon mass ratio for all the four pentaquark states. Since in general the 
Wilson action is known to have 0(a) discretization errors, we extrapolated linearly 
to a — > 0. The mass ratios at hand show remarkably small scaling violations. 
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Figure 4: Continuum extrapolation of the mass ratios m^/ (mjv + my) for the various 
pentaquark states. The horizontal axis shows the nucleon mass in lattice units. The lines 
are linear fits to the data. 

3. Results and conclusions 

Before we present our final results we discuss how we handled the possible sources of 
errors. In lattice QCD there are five typical sources of uncertainties: i. statistical; 
ii. finite size, Hi. quark mass, iv. finite lattice spacing and v. quenching errors. We 
had all of these five error sources under control in our analysis. 
ad i. Statistical uncertainties are well understood and they were included by standard 
techniques. These techniques need uncorrelated configurations, which we ensured by 
separating the analyzed configurations by as much as 1000 sweeps. Statistical errors 
were always estimated using a full jackknife analysis. The statistical uncertainties of 
our pentaquark masses are typically on the few percent level. 

ad ii. Finite size effects in our quenched analysis turned out to be negligible. Our 
spatial box size was around 2 fm. In order to look for finite size corrections we 
also performed calculations on our two coarsest lattices using different box sizes. 
With our statistical uncertainties of order a few percent we could not find any size 
dependence in the hadron masses. 

ad Hi. Quark masses in lattice spectroscopy are usually treated differently in the 
strange and in the light (up or down) sectors. It is straightforward to take the 
physical strange quark mass. Due to computational/technical difficulties most of 
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the analyses use extrapolations from light quark masses, which are larger than the 
physical one. We also applied this method. We included all the extrapolation errors 
into our analysis. 

ad iv. Finite lattice spacing is an inherent source of uncertainty in any lattice analysis. 
We used the standard approach. We performed our calculations at three different 
lattice spacings and extrapolated to the continuum limit. Luckily enough, the mass 
ratios m 5g /(m N + m K ) have a rather weak dependence on the lattice spacing (c.f. 
Fig. |]). Nevertheless, the quite large, order 10% errors on the finally quoted values 
are mostly due to the chiral and the continuum extrapolation. 
ad. v. We used quenched gauge configurations, which is another source of systematic 
errors. Note, however, that in the case of stable hadrons, this is not expected to 



be very important. Indeed, as can be seen from Refs. with an appropriate 
definition of the scale, the mass ratios of stable hadrons are described correctly by 
the quenched approximation on the 1-2% level. We included a conservative estimate 
of 3% quenching uncertainty in our analysis. 

In addition to these five typical sources of uncertainty, there is a somewhat special 
one, connected to the fact that the + mass is rather close to the N+K threshold. 
In the negative parity channel the S-wave scattering state might contaminate the 
spectrum. Though we showed how these states can be successfully separated, more 
work is needed to clarify this issue and finalise the conclusion on the parity. 

Given the fact that in the m -/(m K + m N ) mass ratio we do not observe any 
scaling violations, we could also quote the value of this quantity on our finest lattice, 
which has the smallest error. This would give 

1.073(34), (3.1) 



m K + m N 

or using the physical kaon and nucleon mass 

too- = 1539 ± 50MeV. (3.2) 

Note, that the continuum extrapolation increases the error on the mass ratio, and 
consequently on the mass, by a factor of 3. 

Our final result after continuum extrapolation is shown in Fig. [|. Of the four I p 
channels, the closest in mass to the experimentally observed + is the 0~ state. The 
1~ state can be seen to be about 15% heavier, but still within one standard deviation 
of the + mass. This state, however, which is a member of an isospin triplet, is not 
the + , since the latter is experimentally known to be an isospin singlet |§. The 
parity partner of the 0~ state is almost two times heavier and lies several standard 
deviations above the experimentally observed mass of the + . This suggests that 
the 0~ state should be identified with the + . 

Obviously in this exploratory study we could not perform a high statistics anal- 
ysis. Quenching effects are quite small, nevertheless, it would be advisable to carry 
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Figure 5: The continuum extrapolated mass ratios m^/imN + itik) for the lowest mass 
pentaquark states in the various I p channels. The horizontal line shows the experimentally 
known mass of the + . 

out the calculations with unquenched configurations. We did not make a system- 
atic study of the possible interpolating operators that are likely to have a good 
overlap with the + , either. This study would be also important. Comparing dif- 
ferent spatial and flavour trial wave functions is practically the only way to obtain 
more information about the flavour structure of the + wave function and decide 
among the several possibilities that have already been put forward in the recent 
literature (e.g. one could analyze diquark-diquark-antiquark pictures @). Further- 
more, one expects additional exotic pentaquark hadrons with different spin, isospin 
and strangeness quantum numbers. Lattice QCD might be useful to test these light 
exotic particle cascades numerically. 

In conclusion, we studied spin 1/2 isoscalar and isovector candidates in both par- 
ity channels for the recently discovered + (154O) pentaquark particle in quenched 
lattice QCD. Our analysis took into account all possible uncertainties. The lowest 
mass that we find in the I p = 0~ channel is in complete agreement with the ex- 
perimental value of the + mass. On the other hand, the lowest mass state in the 
opposite parity I p = + channel is much higher. Since the 1 = 1 channel has already 
been excluded by experiment as a candidate for the + ||, this leaves the I p = 0~ 
state to be identified with the + . Our conclusion is consistent with that of the in- 



dependent analysis of S. Sasaki [24[], who separated in addition also the K+N S-wave 
scattering state by observing a second plateau in the effective mass. This successful 
separation is extremely important for establishing the parity of the B + . Clearly, the 
more sophisticated, CPU expensive, cross correlator technique is needed in order to 
unambigously separate the + and the K+N S-wave and its excited states, which 
will finalize the conclusion on the mass and parity. 
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